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•  A  kinetic  model  for  elementary  reactions  of  H2  electro-oxidation  was  developed. 

•  A  new  H2  reaction  mechanism  was  validated  using  pattern  anode  experiments. 

•  Good  agreements  between  simulation  and  experimental  results  were  obtained. 

•  Effect  of  surface  diffusion  proved  to  be  negligible  under  investigated  conditions. 

•  Rate  limiting  steps  were  elucidated,  depending  on  operating  conditions. 
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Good  understanding  of  the  H2  electrochemical  reaction  mechanism  helps  optimize  SOFC  anode  design 
and  improve  its  performance.  Controversies  still  exist  regarding  H2  oxidation  mechanism  despite 
extensive  studies  performed.  In  this  work,  a  new  charge  transfer  reactions  mechanism  proposed  by 
Shishkin  and  Ziegler  (2010)  based  on  Density  Functional  Theory  (DFT)  study  was  investigated  through 
kinetic  modelling  and  pattern  anode  experimental  validation.  The  new  charge  transfer  reactions 
mechanism  considers  hydrogen  oxidation  at  the  interface  of  Ni  and  YSZ.  It  involves  a  hydrogen  atom 
reacting  with  the  oxygen  ions  bound  to  both  Ni  and  YSZ  to  produce  hydroxyl  (Charge  transfer  reaction  1), 
the  latter  reacting  with  the  other  hydrogen  atom  to  form  water  (Charge  transfer  reaction  2).  The  pre¬ 
dictive  capability  of  this  reaction  mechanism  to  represent  our  experimental  results  was  evaluated.  The 
simulated  Tafel  plots  were  compared  with  our  experimental  data  for  a  wide  range  of  H2  and  H20  partial 
pressures  and  at  different  temperatures.  Good  agreements  between  simulation  and  experimental  results 
were  obtained.  Charge  transfer  reaction  1  was  found  to  dominate  the  overall  charge  transfer  reaction 
under  cathodic  polarization.  Under  anodic  polarization,  a  change  in  the  dominating  charge  transfer  re¬ 
action  from  charge  transfer  reaction  1  to  charge  transfer  reaction  2  was  found  when  increasing  the  H20 
partial  pressure. 

©  2013  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

A  solid  oxide  fuel  cell  (SOFC)  is  a  device  that  converts  chemical 
energy  directly  into  electrical  energy  with  high  efficiency.  SOFCs 
can  operate  on  a  variety  of  fuels  and  have  significantly  less  envi¬ 
ronmental  impact  than  the  traditional  combustion  systems  [1—3]. 
Knowledge  of  the  actual  electrochemical  oxidation  occurring  at  the 
anode  is  of  great  importance  for  further  improvement  of  the  anode 
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performance.  However,  even  for  the  simplest  fuel,  H2,  the  detailed 
electrochemical  reactions  are  not  yet  fully  understood  [4-6], 

The  electrochemical  reactions  occur  at  the  triple  phase  bound¬ 
ary  (TPB),  where  gas  phase,  anode  and  electrolyte  join  [4,7],  For 
reaction  kinetic  studies,  pattern  anodes  have  the  advantage  of 
forming  a  well-defined  2-D  structure,  thus  avoiding  the  compli¬ 
cated  3-D  structure  effects  on  the  electrochemical  reactions  that 
take  place  in  cermet  anodes.  Moreover,  influence  of  gas  transport  in 
porous  media  is  eliminated  in  pattern  anodes,  and  the  TPB  length  of 
pattern  anodes  can  be  calculated  accurately  [5,8],  Yttria-Stabilized 
Zirconia  (YSZ)  is  currently  the  most  popular  electrolyte  material 
for  SOFC,  and  Ni/YSZ  pattern  anodes  have  been  commonly  used  in 
the  characterization  of  electrochemical  reaction  mechanisms. 
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In  order  to  elucidate  the  hydrogen  electrochemical  reaction 
mechanism  on  the  Ni/YSZ  anodes,  extensive  experimental  and 
modelling  studies  have  been  performed  so  far,  and  several  reaction 
mechanisms  for  hydrogen  electro-oxidation  have  been  proposed. 

Bieberle  and  Gauckler  [9]  assumed  that  the  charge  transfer  re¬ 
action  is  oxygen  spillover  from  YSZ  surface  to  Ni  surface.  In  their 
work,  a  kinetic  model  describing  a  six-step  elementary  reaction 
mechanism  was  developed  and  validated  using  experimental  re¬ 
sults  from  Bieberle  et  al.  [8],  The  rate-limiting  process  was  claimed 
to  be  either  adsorption/desorption  of  hydrogen  or  removal  of  ox¬ 
ygen  ions  from  YSZ.  In  addition,  both  processes  could  be  coupled 
with  charge  transfer  reaction  and  contributed  to  the  rate¬ 
determining  steps.  However,  effects  of  water  on  the  anode  perfor¬ 
mance  were  evaluated  in  a  small  range  of  water  concentration  for 
both  modelling  and  experimental  studies  (H2O:  0.04-2%  at  700  °C). 
Also,  no  studies  were  performed  under  cathodic  conditions. 
Moreover,  their  experimental  results  varied  considerably  from 
samples  to  samples  [4], 

Vogler  et  al.  [4]  investigated  different  charge  transfer  reactions 
by  comparing  simulation  results  with  the  experimental  data  re¬ 
ported  by  Bieberle  et  al.  [9],  Hydrogen  spillover  with  water  for¬ 
mation  on  YSZ  surface  was  claimed  to  be  the  mechanism  which 
correctly  predicts  the  observed  relationship  between  polarization 
resistance  and  H2  and  H2O  partial  pressures.  They  also  claimed  that 
hydrogen  spillover  to  oxide  ions,  water  dissociation  on  YSZ,  water 
desorption  from  YSZ,  and  surface  diffusion  of  hydroxyl  ions 
adsorbed  on  YSZ  might  be  possible  dominating  processes.  Due  to 
lack  of  experimental  data,  the  simulation  was  performed  within  a 
limited  range  of  conditions  (H2O:  0.04—2%). 

Goodwin  et  al.  [10]  assumed  that  hydrogen  spills  over  to  the  YSZ 
surface  to  form  hydroxyl,  which  subsequently  spills  over  back  to  Ni 
surface  with  water  formation  on  Ni  surface.  They  claimed  that  this 
mechanism  provides  a  good  fit  to  the  experimentally  observed  Tafel 
plots  reported  by  Mizusaki  et  al.  [11],  However,  both  experiment 
and  simulation  were  conducted  in  a  small  range  of  gas  composi¬ 
tions  (H2:  1-10%,  H20:  0.4-1.7%). 

Based  on  density  functional  theory  (DFT),  Shishkin  and  Ziegler 
[12]  proposed  that  the  H2  oxidation  occurs  primarily  on  the  oxygen 
ions  at  the  interface  between  Ni  and  ZrC>2  (or  Y2O3).  DFT  study 
revealed  that  oxygen  vacancy  formation  at  the  interface  of  Ni/YSZ 
was  more  favourable  than  at  the  Ni  or  YSZ  surface.  The  easier  va¬ 
cancy  formation  directly  leads  to  more  negative  charges  trans¬ 
ferring  from  YSZ  to  Ni  surface.  Moreover,  it  was  found  that  the 
energy  barrier  for  water  formation  at  the  interface  is  lower  than 
that  on  Ni  or  YSZ  surface.  These  indicate  that  the  pathway  of  charge 
transfer  reactions  occurring  at  the  interface  is  more  feasible  than 
oxygen  and  hydrogen  spillover  pathways. 

The  main  goal  of  this  study  is  to  perform  a  kinetic  study  based 
on  the  charge  transfer  reactions  proposed  by  Shishkin  and  Ziegler 
(2010)  and  evaluate  the  capability  of  this  model  to  represent 
experimentally  observed  behaviour  of  H2  oxidation  on  Ni/YSZ 
pattern  anodes. 

A  kinetic  model  representing  surface  chemical  reactions,  elec¬ 
trochemical  reactions,  and  surface  diffusion  of  the  H2/H2O  system 
on  a  SOFC  pattern  anode  was  developed  and  built  in  gPROMS. 

Similar  to  the  experimental  procedure,  Tafel  plots  were  simu¬ 
lated  under  a  wide  range  of  H2  partial  pressure  (PH2)  (1  x  104  Pa— 
9  x  104  Pa)  and  H2O  partial  pressure  (PH2O)  (3  x  103  Pa— 5  x  104  Pa) 
conditions.  The  simulated  Tafel  plots  upon  variation  of  pH2  were 
compared  with  the  Tafel  plots  obtained  through  experiments  at 
700  °C.  This  temperature  was  chosen  because  previous  modelling 
results  were  also  reported  at  700  °C,  and  most  of  the  parameters  are 
available  in  the  literature  at  this  temperature.  In  addition,  simu¬ 
lated  Tafel  plots  upon  variation  of  PH2O  were  compared  with 
experimental  Tafel  plots  at  550  °C.  As  presented  in  our  previous 


studies  [13],  stable  structure  and  precise  TPB  length  could  be  ob¬ 
tained  when  the  pattern  anode  was  tested  under  high  amount  of 
H20  (>10%)  conditions  at  550  °C,  but  not  at  700  °C. 

At  700  °C,  thermodynamic  and  kinetic  parameters  for  surface 
chemical  reactions  were  taken  from  Vogler  et  al.  [4[.  Activation 
energies  and  changes  in  Gibbs  free  energy  for  charges  transfer  re¬ 
actions  were  taken  from  Shishkin  and  Ziegler  [12],  Pre-exponential 
factors  for  the  charge  transfer  reactions  were  treated  as  free  fit 
parameters.  At  550  °C,  more  parameters  were  treated  as  free  fit 
parameters  (to  be  discussed  later)  due  to  lack  of  data  in  the  liter¬ 
ature  at  this  temperature. 

2.  Methodology 

2.1.  Experimental 

The  electrolyte  used  in  this  study  is  an  8  mol%  Y2O3  stabilized 
Zr02  single  crystal  disc  with  (100)  orientation  (MTI  corporation), 
and  the  thickness  of  the  electrolyte  is  0.5  mm.  Ni/YSZ  pattern  anodes 
with  Ni  thickness  of  1  pm  were  fabricated.  A  photolithographic 
process  with  bi-layer  resist  system  was  used  for  constructing  Ni 
patterns  on  YSZ  surface.  Details  of  the  fabrication  process  can  be 
found  in  a  previous  study  [13],  Platinum  (Pt)  cathode  and  the  Pt 
reference  electrode  with  thickness  of  100  nm  were  prepared  by  DC 
magnetron  sputtering.  The  Pt  cathode  (not  patterned)  was  placed 
symmetrically  to  the  pattern  anode  on  the  opposite  side  of  the  YSZ 
single  crystal.  A  three  electrodes  configuration  was  applied  to 
characterize  the  Ni/YSZ  pattern  anode  electrochemical  behaviour. 
The  detailed  description  of  the  three  electrodes  and  the  electro¬ 
chemical  test  set-up  can  also  be  found  in  our  previous  study  [13], 
The  cell  was  first  pre-treated  in  H2/H2O  environments  to  achieve 
stable  performance.  Then,  the  Tafel  plots  were  measured  using  a 
Solartron  Electrochemical  Interface  1287  in  the  potentiostatic  mode 
under  different  pH2  (3  x  104  Pa  <  pH2  <  9  x  104  Pa)  and  PH2O 
(3  x  103  Pa  <  pH2  <  5  x  104  Pa)  conditions.  The  imposed  over¬ 
potential  for  Tafel  plots  measurements  were  between  -0.4  and 
0.4  V. 

2.2.  Modelling  approach 

In  traditional  fuel  cells  models,  the  electrochemistry  is  described 
by  global  reactions  and  the  current  density  is  calculated  based  on 
the  Butler-Volmer  equation.  This  simplified  approach  cannot  fully 
represent  the  complex  processes  occurring  at  the  triple  phase 
boundary  [14],  In  this  work,  detailed  1-D  models  representing  the 
coupled  behaviour  of  surface  chemistry,  electrochemistry  and 
surface  diffusion  were  developed.  Similar  models  were  also  used  in 
previous  studies  [4,8,10],  The  model  describes  not  only  the  elec¬ 
trochemical  reactions,  but  also  the  non-faradaic  chemical  reactions 
and  surface-diffusive  transport  on  Ni  and  YSZ  surfaces.  The  reaction 
rates  of  these  chemical  and  electrochemical  reactions  were  ob¬ 
tained  through  mass-action  kinetics.  The  current  density  was 
calculated  directly  via  charge  transfer  reactions.  Both  steady-state 
and  dynamic  electrode  behaviour  can  be  simulated  with  these 
models. 

A  Ni/YSZ  pattern  anode  half-cell  configuration,  which  was  used 
in  both  Goodwin  et  al.  [10]  and  Vogler  et  al.’s  [4]  studies,  was 
chosen  for  this  study.  To  mimic  the  experimental  conditions,  the 
cell  consists  of  a  dense  YSZ  disk  as  the  supporting  electrolyte,  an 
array  of  thin,  straight  Ni  strips  on  one  side  of  the  electrolyte,  a  Pt 
counter  electrode,  and  a  Pt  reference  electrode  on  the  other  side. 
Due  to  the  symmetry  of  the  pattern  anode,  only  one  domain  is 
considered  for  the  simulation.  This  domain  extends  from  the  centre 
line  of  the  Ni  strip  to  the  centre  line  of  the  adjacent  YSZ  strip.  As 
shown  in  Fig.  1,  dNi,  dysz  are  the  width  of  Ni  and  YSZ  strips, 
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Fig.  1.  (a)  Side  view  of  the  pattern  anode,  (b)  the  half-cell  model. 

respectively,  and  fiNi  is  the  thickness  of  the  Ni  layer.  The  red  line  in 
Fig.  1(b)  is  the  TPB  line. 


-  ^(t) 

(5) 

1  =  E  im(hi  -  TSi)  ~  E  -  TSi) 

ieR f,m  ieRr.m 

(6) 

ho.i  +  Cp(T-Tref) 

(7) 

-  a'+|(r+Tre f) 

(8) 

So,i  +  a  (in T  -  lnTref)  +b'(j-  Tref) 

(9) 

Here,  rf,m  and  rr,m  represent  the  forward  and  reverse  reaction 
rates  (mol  s-1  cm-2),  c,-  is  the  area-specific  concentration  of  species 
z  (mol  cm-2),  v'v  v"  are  stoichiometric  coefficients  of  species  i  in  the 
forward  and  reverse  reactions,  k^m  is  the  pre-exponential  factor  of 
the  forward  reaction,  £j?^  is  the  activation  energy  for  the  forward 
reaction  (J  mol-1),  and  /?m  is  a  temperature  coefficient.  AGm  is  the 
Gibbs  free  energy  change  of  the  reaction  (J  mol-1);  hi  (J  mol-1)  and 
Si  0  K-1  mol-1)  are  molar  enthalpy  and  molar  entropy  of  species  i  at 
temperature  T  (K).  /lo.i  and  So.i  are  molar  enthalpy  and  molar  en¬ 
tropy  of  species  i  at  reference  state  temperature  (Tre f :  25  °C),  Cp  is 
the  average  specific  heat  over  the  temperature  range  involved,  a 
and  b  are  constants  used  for  average  specific  heat  calculation. 

The  concentration  change  of  surface  species  i  due  to  chemical 
reactions  is  given  as  [15]: 

(10) 


2.3.  Governing  equations 

2.3.1  Elementary  chemical  reactions 

Both  Ni  and  YSZ  surface  chemical  reactions  in  this  study  are 
written  as  reversible  elementary  reaction  steps.  The  general  form  of 
the  chemical  reaction  is  expressed  as: 

fef,m 

Aj+Bj  ^  Q  +  Dj  (1) 

kr,m 

where  m  represents  all  the  reactions  that  include  adsorption/ 
desorption  of  gas  phase  species,  and  surface  reactions,  A,  Bi,  Q  and 
D,-  are  species  involved  in  these  reactions,  kf,m,  and  krm  are  the 
forward  and  reverse  reaction-rate  constants,  respectively,  and  the 
units  are  given  in  terms  of  mol,  cm  and  s. 

The  reaction  rate  for  forward  and  reverse  reactions  can  be 
represented  using  Equations  (2)— (9).  The  forward  reaction  rate 
constants  are  calculated  using  Arrhenius  expression,  and  the 
reverse  reaction  rate  constants  are  obtained  based  on  a  thermo¬ 
dynamically  consistent  way  [15,16], 

Ff,m  =  kf,m  II  Cf  (2) 


For  surface  chemistry,  the  reaction  rate  is  usually  expressed  in 
the  form  of  surface  coverage  (#;),  which  is  defined  as  [8[: 

6i  =  Ci/Nk  (11) 

»  =  ! 

where  Nk  is  the  area-specific  density  of  Ni  or  YSZ  (mol  cm-2).  The 
surface  coverage  of  each  species  i  is  between  0  and  1,  and  the  total 
surface  coverage  of  all  the  species  is  1.  The  unfilled  surface  sites  are 
also  considered  as  one  individual  species  [10], 

2.3.2.  Charge  transfer  (CT)  reactions 

Charge  transfer  reaction  rates  directly  determine  the  faradaic 
current  produced.  They  are  expressed  in  the  same  way  as  those  for 
the  chemical  reactions  [4,8,17], 

rCT,k  =  kcT,f  n^- ^CT,!  n  (12) 

ie  i?fk  l'eflr,  k 

kcr,f  =  kcr,f  exp  ^  -  q^Jexp  (13) 


k(,m  =  kf,m7'ftne 


rt) 


pKm  n 1 


where  rcr,k  is  the  net  reaction  rate  of  the  charge  transfer  reaction  k 
(mol  s-1  cm-2),  kcr.f,  kcr.r  are  forward  and  reverse  reaction  rate 
constants  (cm3  mol-1  s-1),  kj-j-f  (cm3  mol-1  s-1)  and  £cr,f  0  mol-1) 
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Table  1 

Elementary  reaction  steps  and  kinetic  parameters  used  in  the  work,  which  are  taken  from  Vogler  et  al.  (2009). 


Reversible  reactions 

Pre-exponential  factor  (kj},) 

Activation  energy  (EJf)  (kj  mor1) 

An 

1 

k, 

H2Ogas  +  (YSZ)^  H2Oysz 

6.6  x  1011  cm3  mor1  s-1 

0 

0.5 

2 

Ofe  +  HjOvsz  4  20HySZ 

1.6  x  1022  cm2  mor1  s-1 

9.6 

0 

3 

Vq(YSZ)  +  Ofez  J  Oq(YSZ)  +  (YSZ) 

1.6  x  1022  cm2  mor1  s-1 

90.9 

0 

4 

1C, 

H2Ogas  +  (Ni)  -  H20Ni 

1.4  x  1010  cm3  mol-1  r1 

0 

0.5 

5 

H2  gas  +  2(Ni)  4  2HNi 

9.8  x  1017  cm5  mor1  s-1 

0 

0.5 

6 

ke 

HNi  +  0Ni  7=7  OHNi  +  (Ni) 

5.0  x  1022  cm2  mol-1  s-1 

97.0 

0 

7 

k7 

H20Ni  +  0Nj  20HNi 

5.4  x  1023  cm2  mol  s 

20.9 

0 

8 

kg' 

OHNi  +  HNi  tt:  H20Ni  +  (Ni) 
k-g 

3.0  x  1020  cm2  mol1  s-1 

43.0 

0 

are  pre-exponential  coefficient  and  activation  energy  of  the  for¬ 
ward  reaction  respectively,  AGct  0  mol-1)  is  the  Gibbs  free  energy 
change  of  the  charge  transfer  reaction,  n  is  the  number  of  electrons 
transferred  in  the  electrochemical  reaction,  a  is  the  transfer  coef¬ 
ficient,  r/  is  the  overpotential,  and  F  is  the  Faraday  constant 
(96,487  C  mol-1). 

The  faradaic  current  density  produced  due  to  a  charge  transfer 
reaction  is  described  as  [18]: 

h  =  nFLipBrcTk,  Jf,t  =  ^  h  (15) 

where  fy  is  the  current  density  (A  cm-2),  Ltpb  is  the  TPB  length 
(cm  cm-2).  The  total  current  density  j'f.t  is  the  sum  of  current 
density  produced  by  all  the  charge  transfer  reactions. 

2.3.3.  Surface  diffusion 

Surface  species  diffusion  is  given  by  Fick’s  law,  as  shown  in 
Equation  (16)  [4,15]. 


Table  2 

Surface  diffusion  parameters  taken  from  Vogler  et  al.  (2009). 


Species 

D?  (cm2  s-1) 

Efa  (kj  mor1) 

HNi 

4.6  x  10-3 

14.4 

0Ni 

6.3  x  10-3 

61.3 

OHNi 

6  x  10-3 

30.5 

H20Ni 

6  x  10-3 

28.0 

0?iz 

5.5  x  10-7 

90.0 

H20ysz 

13  x  10-2 

55.0 

ohysz 

1.3  x  10-2 

55.0 

I,  =  06) 

Di  =  Dfexp^-ff^  (17) 


Fig.  2.  Schematic  illustration  of  hydrogen  oxidation  at  the  interface  of  Ni/YSZ. 
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Table  3 

Thermodynamic  parameters  for  surface  species  at  700  °C  and  25  °C. 


Species 

hi  at  25  °C 
(kj  mol-1) 

Si  at  25  °C 

Q  mol-1  K-1) 

hi  at  700  °C 
(kj  mol-1) 

Si  at  700  °C 
(J  mol-1  K 1) 

H2Ogas 

-241.8 

188.8 

— 217.0b 

232.0b 

H2  gas 

0 

130.7 

20.0b 

165.0b 

H2Oysz 

-291.9 

49.3 

-273.0 

98.0 

Ofez 

-183.5 

0a 

-236.0 

0 

OHysz 

— 270.3a 

62.5a 

-283.0 

67.0 

YSZ 

0 

0 

0 

0 

H2Om 

-302.6 

87.6 

-273.0 

130.0 

Hni 

— 48.3a 

28.3 

-32.0 

41.0 

0Ni 

-237.5 

19.4 

-222.0 

39.0 

OHNi 

-214.6 

75.8 

-193.0 

106.0 

Ni 

0 

0 

0 

0 

Vo 

0 

0 

0 

0 

O* 

— 192.0a 

0 

-236.0 

0 

The  other  parameters  at  25  °C  and  700  °C  are  taken  from  Habibzadeh  [19]  and 
Vogler  et  al.  [4],  respectively. 
a  Free  fit  parameters. 
b  Taken  from  O’Hayre  et  al.  [2], 


Kinetic  and  thermodynamic  parameters  used  for  charge  transfer  reactions. 


700  °C 

550  °C 

k§  (cm3  mol  1  s  ') 

3.5  x  1010 

3.5  x  1010 

(cm3  mol-1  s-1) 

2.0  x  1025 

2.0  x  1025 

ACcntkjmol1) 

6.31 

4.30 

AGct2  (kj  mol1) 

-88.8 

-68.7 

AEcri  (kj  mol-1) 

13.0a 

45.0 

AEcn  (kj  mol-1) 

154.9a 

154.9a 

a  Taken  from  Shishkin  and  Ziegler  (2010),  the  others  are  all  fitted  parameters. 


pH2=3.0x!04  Pa,  pH20=3.0xl03  Pa,  T=700°C 


Potential  vs.  reference  potential  /  V 


pH20=3.0x!03  Pa,  T=700°C 


Potential  vs.  reference  potential  /  V 
Fig.  4.  Effect  of  pH2  on  the  current  density  at  700  °C. 


where  J;  is  the  surface  flux  of  adsorbate  i  (mol  cm-2  s-1),  q  is  the 
concentration  of  species  i  (mol  cm-3),  D,  is  the  surface  diffusion 
coefficient  of  species  i  (cm2  s-1),  z  is  the  distance  along  the  surface 
perpendicular  to  the  TPB  line  (cm),  Df  is  the  pre-exponential 
constant  (cm2  s-1),  and  £frt  is  the  activation  energy  (J  mol-1). 


pH2=5x!04  Pa,  pH20=3  xl03Pa,  T=700°C 


Potential  vs.  reference  potential  /  V 


pH2=7x!04  Pa,  pH20=3  xl03Pa,  T=700°C 


Potential  vs.  reference  potential  /  V 


pH2=9x!04  Pa,  pH20=3  xl03Pa,  T=700°C 


Potential  vs.  reference  potential  /  V 


Fig.  3.  Comparison  between  experimental  and  simulation  Tafel  plots  upon  variation  of  pH2  at  700  °C. 
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pH2=2.0xl04  Pa,  pH20=3.0xl03  Pa,  T=550°C 


Potential  vs.  reference  potential  /  V 


pH2=2.0x!04  Pa,  pH20=1.0x!04  Pa,  T=550°C 


Potential  vs.  reference  potential  /  V 

pH2=2.0x!04  Pa,  pH20=5.0x!04  Pa,  T=550°C 


Potential  vs.  reference  potential  /  V  „  .  ,  .  ,  , „ , 

Potential  vs.  reference  potential  /  V 

Fig.  5.  Comparison  between  experimental  and  simulated  Tafel  plots  upon  variation  of  pH20  at  550  °C. 


The  surface  concentration  change  for  each  adsorbed  species  i, 
which  is  described  in  Equation  (18),  is  determined  by  both  chemical 
reaction  and  surface  diffusion  [4], 


Because  of  the  symmetry  of  the  model,  the  boundary  conditions 
re  [10]: 

At  z  —  0.5di\ii  and  z=  0.5dvsz, 


At  the  TPB  line,  for  reactant  species  directly  involved  ii 
charge  transfer  reactions, 


For  product  species  directly  involved  in  the  charge  transfer 
reactions, 


For  those  surface  species  that  are  not  involved  in  charge  transfer 
reactions, 


2.4.  Investigation  of  hydrogen  oxidation  at  Ni/YSZ  interface 

In  this  study,  chemical  reactions  on  the  Ni  surface  were  referred 
from  Bieberle  and  Gauckler  (2002),  Vogler  et  al.  (2009)  and 
Goodwin  et  al.  (2009).  The  YSZ  surface  reactions  were  the  same  as 


PH2=2.0xl04  Pa,  T=550°C 


Fig.  6.  Effect  of  pH20  on  the  cu 


Ji  =  0 


:  density  at  550  "C 
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pH20=3.0xl03  Pa,  T=550°C  pH2=2.0x!04  Pa,  T=700°C 


Potential  vs.  reference  potential  /  V  Potential  vs.  reference  potential  /  V 


Fig.  7.  Effect  of  pH2  and  pH20  on  the  current  density  at  550  °C  and  700  °C,  respectively. 


used  in  Vogler  et  al.  [4],  Charge  transfer  reactions  are  based  on  DFT 
study  performed  by  Shishkin  and  Ziegler  [12],  Table  1  summarizes 
the  elementary  reaction  steps  on  Ni  and  YSZ  surfaces,  as  well  as  the 
corresponding  kinetic  parameters  (pre-exponential  factors  and 
activation  energy  for  forward  reactions  at  700  °C  used  in  this  work). 

Regarding  the  charge  transfer  reactions,  Shishkin  and  Ziegler 
[12]  claimed  that  oxygen  vacancy  formation  at  the  interface  of  Ni/ 
YSZ  is  more  favourable  than  at  the  YSZ  surface,  which  results  in 
more  negative  charges  being  transferred  from  YSZ  to  Ni  surface. 
Moreover,  it  was  found  that  the  energy  barrier  of  water  formation 
at  the  interface  was  lower  than  that  on  Ni  or  YSZ  surface.  These 
indicate  that  charge  transfer  reactions  occurring  at  the  interface  are 
more  likely  than  oxygen  and  hydrogen  spillover.  In  addition,  water 
produced  at  the  interface  is  assumed  to  be  in  the  gas  phase,  since 
the  energy  barrier  for  water  desorption  from  the  interface  is  very 
small,  based  on  DFT  calculation  [12],  The  charge  transfer  reactions 
are  listed  in  reactions  (23)  and  (24),  and  schematically  demon¬ 
strated  in  Fig.  2.  H2  oxidation  at  the  Ni/YSZ  interface  involves  a 
hydrogen  atom  reacting  with  the  oxygen  ions  bound  to  both  Ni  and 
YSZ  surfaces  to  produce  hydroxyl  ion,  which  reacts  with  the  other 
hydrogen  atom  to  form  water. 

HNi  +  Ofnter  ^  OHinter  +  (Ni)  +  eNi  (23 ) 

HNi  +  OHinter  -  H2Ogas  +  (inter)  +  (Ni)  +  eNi  (24) 

In  this  reaction  mechanism,  there  are  two  gas  phase  species  (H2, 
H2O),  four  Ni  surface  species  (HNi,  Oni,  OHm,  FhONi),  three  YSZ 
surface  species  (F^Oysz,  Oysz,  OHysz),  two  bulk  species  (lattice 
oxygen  (Oq)  and  oxygen  vacancy  (V5)),  and  three  species  at  the 
interface  (0^er,  OHr,ter,  and  free  surface  (inter)).  The  Ni  free  sur¬ 
face  (Ni)  and  YSZ  free  surface  (YSZ)  are  also  considered  as  two 
single  species.  At  700  °C,  Shishkin  and  Ziegler  [12]  calculated  the 
activation  energy  for  the  forward  reactions  of  the  charge  transfer 
reactions,  and  found  13  kj  mol-1  and  154.9  kj  mol-1  for  charge 
reactions  (23)  and  (24),  respectively.  Changes  in  Gibbs  free  energy 
of  the  charge  transfer  reactions  obtained  through  DFT  calculation 
are  113  kj  mol-1  and  -75.4  kj  mol-1,  respectively.  Surface  diffusion 
parameters  are  taken  from  Vogler  et  al.  [4]  and  are  listed  in  Table  2. 

3.  Results  and  discussion 

Thermodynamic  parameters  for  surface  species  at  25  °C  and 
700  °C  used  in  the  model  are  shown  in  Table  3.  Kinetic  and 


thermodynamic  parameters  for  charge  transfer  reactions  were 
either  taken  from  Shishkin  and  Ziegler  or  determined  through 
fitting  the  experimental  results.  Those  parameters  are  given  in 
Table  4. 

At  700  °C,  the  activation  energy  of  charge  transfer  reactions 
(Reactions  (23)  and  (24))  were  taken  from  Shishkin  and  Ziegler 
[12],  It  was  realized  that  the  location  for  the  minimum  of  the  log¬ 
arithm  of  the  current  density  on  the  Tafel  plots  depends  on  the 
Gibbs  free  energy  change  of  both  charge  transfer  reactions.  This 
Gibbs  free  energy  changes  were  then  adjusted  for  best  fit,  resulting 
in  a  value  of  6.31  kj  mol-1  for  Reaction  (23)  and  a  value 
of  -88.8  kj  mol-1  for  Reaction  (24).  It  is  also  important  to  notice 
that  one  would  expect  a  current  density  close  to  zero  at  OCV,  but 
the  experimental  and  simulation  results  give  a  baseline  current.  For 
the  experimental  data,  the  reason  for  a  baseline  current  is  due  to 
the  chemical  potential  difference  between  the  cathode  and  anode, 
which  leads  to  a  small  exchange  current  density  at  equilibrium 
conditions.  Gas  leaks,  as  well  as  electronic  leaks,  may  also  cause  a 
small  exchange  current  density.  For  the  simulation,  the  value  of  the 
minimum  current  density  is  extremely  sensitive  to  very  small 
changes  in  the  value  of  the  Gibbs  free  energy.  For  example,  in 
Fig.  3(b),  the  set  of  Gibbs  free  energies  yields  a  very  small  current 
density  (log(i)  =  -10.76  A  cm-2).  Would  we  change  the  value  of 
those  free  energies  by  less  than  0.5%,  the  logarithm  of  the  minimum 
current  density  would  increase  to  -6.0  A  cm-2.  Note  that  changes 
in  Gibbs  free  energy  of  this  magnitude,  although  they  did  affect 
significantly  the  minimum  value  of  the  current  density,  they  did  not 
have  any  effect  on  the  values  of  the  cathodic  and  anodic  current 
densities.  Also,  note  that  the  value  of  6.31  kj  mol-1  for  Reaction  (23) 
is  significantly  smaller  than  the  calculated  value  through  DFT 
(113  kj  mol-1)  by  Shishkin  and  Ziegler.  The  Gibbs  free  energy 
change  of  the  overall  H2  oxidation  reaction  is  the  sum  of  the  Gibbs 
free  energy  change  of  all  the  elementary  reactions.  With  the  fitting 
value  (6.31  kj  mol-1 ),  the  sum  of  the  Gibbs  free  energy  change  of  all 
the  elementary  reactions  is  much  closer  to  the  Gibbs  free  energy 
change  of  the  overall  H2  reaction  than  when  using  the  DFT  calcu¬ 
lated  value. 

At  550  °C,  charges  in  Gibbs  free  energy  as  well  as  the  activation 
energy  for  Reaction  (23)  were  treated  as  a  free  fit  parameter  to 
obtain  good  fitting  between  simulation  and  experimental  results. 
The  activation  energy  for  Reaction  (24)  had  insignificant  effect  on 
the  simulation  results,  and  was  taken  from  the  DFT  calculation. 

In  Fig.  3,  both  simulation  and  experimental  results 
regarding  the  effect  of  H2  partial  pressure  on  current  density 
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are  presented.  The  simulation  and  experimental  conditions  are: 
3  x  104  Pa  <  pH2  <  9  x  104  Pa,  pH20  =  3  x  103  Pa, 
T=700°C. 

In  Fig.  3  (as  well  as  in  Fig.  5)  the  simulation  results  from  models 
with  and  without  surface  diffusion  overlap  exactly,  which  means 
that  the  surface  diffusion  is  fast  under  the  above  operating  condi¬ 
tions.  On  both  Ni  and  YSZ  surfaces,  surface  coverage  profiles  are 
uniform.  The  experimental  results  can  be  well  reproduced  by 
simulation  under  each  pH2  condition.  Fig.  4  shows  that  pH2  has 
positive  effects  on  the  current  density  under  anodic  polarization, 
and  the  effect  of  pH2  on  the  cathodic  current  density  is  less  pro¬ 
nounced,  which  agrees  well  with  the  experimental  data.  As  seen  in 
Fig.  4,  the  experimentally  observed  anodic  current  density  in¬ 
creases  when  increasing  pH2,  and  the  experimental  cathodic  cur¬ 
rent  density  is  almost  independent  of  pH2. 

Fig.  5  evaluates  the  influence  of  pH20  on  the  current  density  for 
both  experiments  and  simulation.  The  simulation  and  experimental 
conditions  are:  3  x  103  Pa  <  pH20  <  5  x  104  Pa,  pH2  =  2  x  104  Pa, 
T  =  550  °C.  The  effect  of  surface  diffusion  on  both  anodic  and 
cathodic  current  density  is  negligible.  Under  high  cathodic  poten¬ 
tial,  the  simulated  current  density  is  larger  than  the  experimental 
current  density,  but  simulation  results  still  reasonably  predict  the 
experimental  data.  As  shown  in  Fig.  6,  the  experimentally  observed 
decrease  in  cathodic  current  density  with  increasing  pH20  can  be 
reproduced  by  the  simulation.  For  both  modelling  and  experi¬ 
mental  results,  the  effect  of  pH20  on  anodic  current  density  is 
small. 

In  order  to  compare  the  simulation  results  with  the  experi¬ 
mental  results,  the  effect  of  pH2  on  current  densities  was  simulated 
at  700  °C  and  the  effect  of  pH20  on  current  densities  was  simulated 
at  550  °C,  which  are  presented  in  Fig.  7.  As  mentioned  earlier,  it  was 
not  possible  to  conduct  a  complete  kinetic  study  experimentally  at 
700  °C  because  the  Ni/YSZ  pattern  anodes  experienced  significant 
structure  changes  at  this  temperature  for  H20  greater  than  10%. 
Nonetheless,  effect  of  H20  at  700  °C  can  be  studied  from  simulation. 
Fig.  7  shows  the  dependence  of  simulated  current  density  on  pH20 
at  700  °C  as  well  as  the  effect  of  pH2  on  current  density  at  550  °C. 
The  effect  of  pH2  on  current  density  at  550  °C  could  have  been  done 
experimentally,  but  it  was  deemed  sufficient  here  to  conduct  this 
study  of  pH2  effect  on  current  density  only  with  simulated  results. 
At  both  550  °C  and  700  °C,  surface  diffusion  does  not  show  any 
effects  on  the  current  density  under  different  pH2  and  pH20  con¬ 
ditions.  At  550  °C,  pH2  has  a  positive  effect  on  the  anodic  current 
density  and  a  very  weak  effect  on  cathodic  current  density,  which  is 
similar  to  the  simulation  results  obtained  at  700  °C.  At  700  °C,  the 
cathodic  current  density  increases  when  increasing  pH20,  whereas 
the  dependence  of  anodic  current  density  on  pH20  is  much  weaker. 
This  is  also  consistent  with  the  simulation  results  obtained  at 
550  °C. 

Table  5  presents  the  surface  coverage  of  Ni  and  YSZ  surface 
species  for  various  pH2  and  pH20  conditions  when  an  overpotential 
of  0.3  V  is  imposed.  When  increasing  both  pH2  and  pH20,  Ni  free 
sites  slightly  decrease.  With  an  increase  in  pH2,  #(HNi)  increases,  but 
0(ONi)  and  0(OHNi)  decrease.  Surface  coverage  of  H2Oni  and  YSZ 
surface  species  are  slightly  affected  by  pH2.  0(Oysz)  slightly  de¬ 
creases,  and  0(OHySZ)  increases.  When  increasing  pH20,  surface 
coverage  of  Oysz  decreases  and  more  H2Oysz  and  OHySZ  are 
generated.  In  addition,  the  amount  of  H2C>Ni,  Om  and  OH^  increase, 
and  Hnj  almost  remains  constant  with  increasing  pH20.  At  the  Ni/ 
YSZ  interface,  surface  coverage  of  OHj^ter  slightly  increases  and 
^(Ojnter)  shows  small  decrease  with  increasing  both  pH2  and  pH20. 

As  discussed  before,  the  anodic  current  density  increases  with 
an  increase  in  pH2.  Therefore,  the  trend  of  current  density  change  is 
related  to  the  change  in  0(HNi)  and  0(OH^ter)  upon  variation  of  pH2 
under  constant  pH20  and  constant  temperature.  When  increasing 
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pH2=5.0xl04  Pa,  pH20=3.0xl03  Pa,  T=700°C 


Potential  vs.  reference  potential  /  V  Potential  vs.  reference  potential  /  V 


Potential  vs.  reference  potential  /  V 
Fig.  8.  Dependence  of  surface  coverages  of  YSZ  surface  species  and  current  i 


Potential  vs.  reference  potential  /  V 
i  electrode  potential:  pH2  =  5.0  x  104  Pa,  pH20  =  3.0  x  103  Pa,  T  =  700  °C. 


pH2,  0(Hni)  and  0(OH7,ter)  increases  and  leads  to  an  increase  in 
anodic  current  density.  Under  cathodic  polarization,  the  extent  of 
0(HNj)  change  is  the  same  as  that  under  anodic  polarization  when 
increasing  the  same  amount  of  H2.  But,  the  extent  of  0(OFIj^ter) 
change  under  cathodic  polarization  is  smaller  than  under  anodic 
polarization.  This  explains  the  strong  dependence  of  the  anodic 
current  density  and  the  weak  dependence  of  the  cathodic  current 
density  on  pH2.  When  pH2  and  temperature  are  kept  constant,  the 
current  density  increases  when  increasing  pH20,  and  the  trend  of 
current  density  changes  is  mainly  related  to  ^(OHinter)  changes. 
According  to  the  charge  transfer  reactions  (23)  and  (24),  electrons 
produced  are  directly  related  to  0(HnO,  0(OH£,ter)  and  0(Ofn^er). 
Under  constant  pH2  and  constant  temperature,  0(HNi)  is  almost 
constant,  and  0(O|^er)  decreases  with  an  increase  in  pH20  leading 
to  more  0(OHj^ter)  production,  which  explains  the  increase  in  cur¬ 
rent  density.  In  addition,  the  extent  of  @(OHinter)  change  is  bigger 
under  cathodic  polarization  than  under  anodic  polarization  when 
increasing  the  same  amount  of  H20.  This  is  the  reason  why  the 
cathodic  current  density  depends  more  strongly  on  pH20  than  the 
anodic  current  density  does. 

Under  constant  pH2  and  pH20  conditions,  0(H2Oni)  and  0(HNi) 
increase  with  decreasing  temperature,  while  surface  coverages  of 
other  Ni  surface  species  decrease.  On  the  YSZ  surface,  0(H2Oysz)  and 
0(OHysz)  increases  when  decreasing  the  temperature.  Surface 
coverages  of  other  YSZ  surface  species  remain  constant  with  a 
decrease  in  temperature.  At  the  interface,  0(OHr,ter)  increases  when 
decreasing  temperature. 


On  Ni  surface,  surface  coverage  of  different  species  are  not  influ¬ 
enced  by  varying  the  electrode  potential.  But  surface  coverage  of 
species  on  YSZ  surface  slightly  change  when  changing  the  electrode 
potential.  As  shown  in  Fig.  8(a)-(c),  from  cathodic  polarization  to 
anodic  polarization,  surface  coverage  of  OHySz  increases  by  about 
1.78%,  and  the  surface  coverage  of  other  YSZ  surface  species  slightly 
decrease.  ^(Oysz)  decreases  by  about  0.0013%,  and  surface  coverage  of 
free  YSZ  surface  decreases  by  about  0.005%.  Although  the  changes  of 
6(0 ysz)  and  0(YSZ)  are  very  small  and  almost  negligible,  a  noticeable 
difference  in  current  density  is  caused.  Fig.  8(d)  describes  the  effect  of 
electrode  potential  on  total  current  density  ( j )  as  well  as  current 
density  contributed  by  charge  transfer  reaction  (23)  (j i)  and  (24)  (/2). 
The  total  current  density  is  the  sum  of  jj  and  j2.  Fig.  8(d)  shows  that  the 
total  current  density  is  mainly  due  to  ji  under  the  operating  condi¬ 
tions  for  this  figure.  The  contribution  of  current  density  from  reaction 
(24)  is  small.  Therefore,  the  charge  transfer  reaction  (23)  dominates 
the  overall  charge  transfer  reaction  for  H2  oxidation. 

Under  anodic  polarization,  the  trend  of  0(OHySz)  change  with 
electrode  potential  is  consistent  with  the  trend  of  current  density 
changes  when  pH2  and  pH20  are  kept  constant.  The  overall  H2 
oxidation  becomes  faster  with  an  increase  in  0(OHySZ). 

Under  cathodic  polarization,  current  density  change  with  elec¬ 
trode  potential  follow  the  trend  of  changes  in  0(Oysz)  and  free  YSZ 
surface.  With  more  free  YSZ  surface,  the  overall  H2  oxidation  ac¬ 
celerates  and  the  generated  current  density  increases.  Similar 
behaviour  was  observed  when  changing  pH2,  while  keeping  pH20 
(3.0  x  103  Pa)  and  temperature  (700  °C)  constant. 
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Fig.  9.  Dependence  of  surface  coverages  of  YSZ  surface  species  and  current  i 
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;  on  electrode  potential:  pH2  =  2.0  x  104  Pa,  pH20  =  2.0  x  104  Pa,  T=  550  °C. 


Fig.  9  shows  the  changes  at  550  °C  in  various  surface  species  and 
current  density  when  varying  the  electrode  potential.  At  550  °C,  the 
effect  of  electrode  potential  on  surface  coverages  of  YSZ  surface 
species  and  current  density  was  evaluated  under  different  pH20 
conditions.  When  pH20  is  3.0  x  103  Pa,  the  dependence  of  surface 
coverages  and  current  density  on  electrode  potential  is  very  similar 
to  what  was  observed  at  700  °C  under  the  same  pH20  condition 
(not  shown  here).  With  further  increase  in  pH20  to  2.0  x  104  Pa, 
surface  coverage  of  OHysz  increases  and  surface  coverages  of 
H2Oysz  and  free  YSZ  surface  decrease  with  a  decrease  in  cathodic 
potential.  Surface  coverages  of  all  the  species  remain  constant  un¬ 
der  anodic  polarization  (see  Fig.  9(a)— (c)).  In  Fig.  9(d),  the  total 
cathodic  current  density  is  mainly  due  to  ji,  and  j2  has  minor  effect 
on  the  total  cathodic  current  density.  Under  anodic  polarization, 
is  similar  to  j2  and  the  current  density  from  both  charge  transfer 
reactions  equally  affects  the  total  current  density.  Therefore,  charge 
transfer  reaction  (23)  dominates  the  whole  charge  transfer  reaction 
under  cathodic  polarization,  and  both  charge  transfer  reactions 
((23)  and  (24))  are  rate  controlling  under  anodic  polarization. 

In  Fig.  10,  changes  in  surface  coverage  as  well  as  current  density 
upon  variation  of  electrode  potential  are  presented  when  the  par¬ 
tial  pressure  of  H20  is  5.0  x  104  Pa.  Under  cathodic  conditions, 
surface  coverage  of  OHysz  increases  and  surface  coverages  of  Oysz. 
H2Oysz  and  free  YSZ  surface  decrease  with  decreasing  cathodic 
potential,  which  is  similar  to  surface  coverage  changes  under  other 
pH20  conditions.  The  current  density  from  charge  transfer  reaction 
(23)  is  still  the  main  contributor  to  the  total  current  density.  Under 
anodic  polarization,  surface  coverage  of  OH2sz  decreases  and 


surface  coverages  of  Oysz.  H2Oysz  and  free  YSZ  surface  increase 
when  decreasing  anodic  potential,  which  is  opposite  to  surface 
coverage  changes  when  pH20  is  3.0  x  103  Pa.  The  total  anodic 
current  density  is  mainly  from  charge  transfer  reaction  (24),  and 
charge  transfer  reaction  (23)  has  minor  influence  on  the  total 
anodic  current  density.  Charge  transfer  reaction  (24)  dominates  the 
overall  charge  transfer  reaction. 

Under  cathodic  polarization,  the  charge  transfer  reaction  (23) 
proceeds  towards  generating  FI  and  02~  species,  which  leads  to 
an  increase  in  d(02~)  and  a  decrease  in  0(OFT)  when  increasing  the 
cathodic  potential. 

Under  anodic  polarization,  the  charge  transfer  reaction  (23) 
proceeds  towards  the  production  of  hydroxyl  ions  and  electrons, 
which  results  in  an  increase  in  0(OH“)  and  a  decrease  in  0(O2-) 
when  decreasing  the  anodic  potential  at  low  value  of  pH20  is  low 
(3.0  x  103  Pa).  When  increasing  pH20,  more  hydroxyl  ions  are 
produced  through  reaction  2,  as  seen  in  Table  1,  which  increases  the 
reaction  rate  of  charge  transfer  reaction  (24).  The  increase  of  OH- 
also  leads  to  more  O2  production  through  reaction  (23),  which 
explains  the  increase  of  0(O2-)  under  anodic  polarization  when 
pH20  is  high  (5.0  x  104  Pa). 

At  700  °C,  simulated  results  similar  to  what  was  observed  at 
550  °C  were  obtained  when  varying  pH20. 

The  good  agreement  between  the  simulation  and  experimental 
data  in  a  wide  range  of  pFI2  and  pH20  conditions  reveals  that  the 
new  charge  transfer  reactions  mechanism  proposed  by  Shishkin 
and  Ziegler  [12]  can  reasonably  predict  the  experimentally 
observed  behaviour  of  Ni/YSZ  pattern  anode  in  H2/H20 
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pH2=2.0xl04  Pa,  pH20=5.0x!04  Pa,  T=550°C 


Potential  vs.  reference  potential  /  V 


Potential  vs.  reference  potential  /  V 


Potential  vs.  reference  potential  /  V  Potential  vs.  reference  potential  /  V 

Fig.  10.  Dependence  of  surface  coverages  of  YSZ  surface  species  and  current  densities  on  electrode  potential:  pH2  =  2.0  x  104  Pa,  pH20  =  5.0  x  104  Pa,  T  =  550  °C. 


environments.  Surface  diffusion  is  fast  and  does  not  affect  the 
pattern  anode  performance.  Under  cathodic  polarization,  the 
dominating  charge  transfer  reaction  (reaction  (23))  does  not 
change  when  varying  the  temperature  and  partial  pressures  of  H2 
and  H2O.  Under  anodic  polarization,  the  dominating  charge  trans¬ 
fer  reaction  (still  reaction  (23))  does  not  change  when  changing 
pH2  and  temperature.  However,  it  changes  from  charge  transfer 
reaction  (23)  to  charge  transfer  reaction  (24)  when  increasing  the 
partial  pressure  of  H20  from  3.0  x  103  Pa  to  5.0  x  104  Pa. 


4.  Conclusions 


Experimental  data  (both  anodic  and  cathodic  currents  in  the 
Tafel  plots)  can  be  well  reproduced  by  simulation  under  a  wide 
range  of  pH2  and  PH2O  conditions. 

It  was  found  that  the  effect  of  surface  diffusion  on  the  pattern 
anode  performance  is  negligible  under  all  conditions  in  this  work. 
Under  cathodic  polarization,  charge  transfer  reaction  1  is  always 
the  dominating  charge  transfer  reaction,  irrespective  of  pH2, 
PH2O  and  temperature. 

Under  anodic  polarization,  the  dominating  charge  transfer  re¬ 
action  switches  from  charge  transfer  reaction  1  to  charge 
transfer  reaction  2  when  increasing  the  partial  pressure  of  H2O 
from  3.0  x  103  Pa  to  5.0  x  104  Pa. 


This  work  presents  a  kinetic  modelling  study  of  H2  electro¬ 
chemical  oxidation  using  Ni/YSZ  pattern  anodes  as  model  anodes. 
The  model  represents  elementary  reactions  of  H2  oxidation  on  Ni/ 
YSZ  pattern  anodes.  Surface  reactions,  electrochemical  reactions  and 
surface  diffusions  are  described  in  the  model.  The  charge  transfer 
reactions  mechanism  proposed  by  Shishkin  and  Ziegler  [12]  was 
implemented  into  the  kinetic  model  and  applied  to  interpret  our 
experimental  results.  The  two  charge  transfer  reactions  are: 

Charge  transfer  reaction  1  :  HNi  +  ofnter^OHjnter  +  (Ni)  +  eNi 
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Nomenclature 


Charge  transfer  reaction  2  :  HNi  +  OHinter^H2Ogas  +  (inter) 
+  (Ni)  +  eNi 

The  main  results  are  summarized  as  follows. 


Roman  symbols 

At,  Bj,  Q,  D;  species  involved  chemical  reactions 
a',  b'  constants 
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q  area-specific  concentration  of  surface  species  i 

(mol  cm-2) 

D,-  diffusion  coefficient  of  surface  species  i  (cm2  s-1) 

Dp  pre-exponential  constants  (cm2  s  _1) 
deiyt  electrolyte  thickness  (cm) 

Ecr.f  activation  energy  for  forward  charge  transfer 

reaction  (kj  mol-1) 

Efct  activation  energy  for  surface  diffusion  (kj  mol-1) 

E  Faraday  constant  (96487  C  mol-1 ) 

AG  Gibbs  free  energy  change  of  the  reaction  (kj  mol-1) 
hi  molar  enthalpy  (kj  mol-1) 

ho,,  molar  enthalpy  at  reference  state  temperature  (kj  mol-1) 
j  current  density  (A  cm-2) 

jf  Faradic  current  density  (A  cm-2) 

j'pt  total  faradic  current  density  (A  cm-2) 

Ji  molar  flux  of  surface  species  i  (mol  cm-1  s-1) 

kf,m  forward  reaction  rate  constant  (mol,  cm,  s) 
kr,m  reverse  reaction  rate  constant  (mol,  cm,  s) 
kjm  pre-exponential  coefficient  for  the  forward 

reaction  (mol,  cm,  s) 

n  electron  number  transferred  during  an  electrochemical 

reaction 

Nk  area-specific  density  (mol  cm-2) 

p  pressure  (Pa) 

r  reaction  rate  (mol  cm-2  s-1) 

R  gas  constant  (8.314  J  K-1  mol-1) 

Rp  polarization  resistance  (Q  cm2) 

Si  molar  entropy  0  mol-1  K-1) 

So  molar  entropy  at  reference  state  temperature 

0  mol-1  K-1) 

T  temperature  (K) 

v\  stoichiometric  factor  of  species  i  in  forward  reaction 

v'l  stoichiometric  factor  of  species  i  in  reverse  reaction 

z  distance  along  the  surface  perpendicular  to  the  TPB 

line  (cm) 


Greek  symbols 

a  transfer  coefficient 

t]  overpotential  (V) 

6  surface  coverage 


Subscripts 

act  activation 


CT  charge  transfer  reaction 

CTRs  charge  transfer  reactions 

elyt  electrolyte 

F  Faraday 

f,m  forward  reaction 

H2  hydrogen 

i  surface  species 

k  Ni  or  YSZ  surface 

m  adsorption,  desorption  and  surface  reactions 

r,m  reverse  reaction 


Abbreviation 

DFT  density  functional  theory 

Ni  nickel 

SOFC  solid  oxide  fuel  cell 

TPB  triple  phase  boundary 

YSZ  yttria-stabilised-zirconia 
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